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We discuss the thermodynamics of two-colour QCD with four flavours of staggered quarks on S'' x 4 and 16"' x 4 
lattices. In our simulations we use the Naik action for the fermions and a (1,2) tree-level improved gauge action. We 
analyze the deconfinement and chiral phase transitions for four different quark masses (m=0. 1,0.05,0.025,0.015). 
Contrary to three-colour QCD the peak in the Polyakov loop susceptibility decreases with decreasing quark mass. 
This reflects an early breaking of the string in the heavy quark potential, which we verify explicitly by calculating 
the heavy quark potential at finite temperature using Polyakov loop correlations. 



1. Introduction 

We investigate the thermodynamics of QCD 
with four flavours and colour group SU (2) in the 
staggered discretisation [|l|. In the continuum 
limit the phase transition is expected to be of 
first order induced by fluctuations 

In our simulations we used the exact hybrid- 
MC algorithm with a Symanzik improved (1,2) 
gluon action and the Naik improvement for the 
fermionic sector. To reduce the number of conju- 
gate gradient steps we used the minimal residual 
extrapolation method (MRE)|^. The residuum 
for the CG was 10^'* inside the molecular dy- 
namic trajectory and 10^^^ elsewhere. Each tra- 
jectory was of length 1, while the number of steps 
in the molecular dynamic evolution was adjusted 
to get a good acceptance. 

For the analysis of the phase transition lat- 
tices of size 8'^ X 4 with three different bare quark 
masses m — 0.1, 0.05, 0.025 and size 16^ x 4 with 
m = 0.015 were used to calculate the Polyakov 
loop, the chiral condensate and their susceptibil- 
ities. We used a noisy estimator with 25 random 
sources to calculate the susceptibility of V'V'- Sim- 
ulations on 16^ X 4 lattices at m = 0.025 were 
performed to calculate finite temperature poten- 
tials from Polyakov loop correlations. For com- 
parison of these potentials we also used results of 
two flavour QCD with colour group SU{3) and 
standard gluonic and fermionic actions B . 



2. Polyakov Loop 

The results for the Polyakov loop are shown 
in Fig.l. With decreasing mass they seem to get 
flatter showing no clear sign of a phase-transition. 
For small /3-values the Polyakov loop does not 
drop to zero which is a first hint for string break- 
ing. Comparing the data for the susceptibility of 
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Figure 1. Polyakov loop on 8"^ x 4 lattices (m=0.1, 
0.05, 0.025) and on 16^ x 4 lattices (m=0.015). 

the Polyakov loop for SU{2) with four fiavours 
(Fig. 2) and SU{3) with two flavours [|| one sees 
quite a different mass dependence of the peak 
heights. For SU{2) the peaks decrease with de- 
creasing mass while for SU{3) the heights have 
the tendency to increase slightly. This reflects a 
smooth crossover from the low to the high tem- 
perature phase and an early breaking of the string 
for 5*1/(2) (see Fig. 5). Further simulations are 
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Figure 2. Susceptibility of the Polyakov loop on 
8^ X 4 lattices (m=0.1, 0.05, 0.025) and on 16^ x 4 
lattices (m=0.015). 

needed to see if this behaviour is also valid for 
larger Nj- or if it could be a strong coupling ef- 
fect. 

3. Chiral Condensate 
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Figure 3. Chiral condensate for m=0.1, 0.05, 
0.025 (8^ X 4), m=0.015 (16^ x 4) and the lin- 
ear extrapolation to m=0. 

The data for the chiral condensate in Fig. 3 
show a clear phase transition in the chiral limit 
although (ipip) stays finite for finite bare quark 
masses. The results for the two smallest masses 
were used to extrapolate to the zero mass limit 
with a linear ansatz 

{ipip)(m = 0) = {'>p'ip){m) — c ■ m. (1) 
For P > 1.2 chiral symmetry is restored in the 



zero mass limit within errors. A wider range of 
masses are needed to use other ansatze for the 
extrapolation. 

So far, our results are compatible with earlier 
investigations by Kogut et al. 0. In contrast to 
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Figure 4. Susceptibility of the chiral condensate 
form=0.1, 0.05, 0.025 (8^x4), m=0.015 (16^x4). 

the susceptibility of L the comparison of SU{2) 
and SU{S) data |^ shows no major difference for 
the susceptibility of the chiral condensate. Both 
exhibit clear peaks that increase with decreasing 
bare quark mass. Also quantitatively the quark 
mass dependence of the peak heights is quite sim- 
ilar for both gauge groups. We have fitted the 
peak heights to the ansatz 



(2) 



Action 


Lattice 




SU{2),Nf = 4 
SU{3),Nf = 2 
SU{3),Nf = 2 


8^ X 4 
8^ X 4 
16^ X 4 


0.95(7) 
0.84(5) 
0.93(8) 



Table 1. Fit-results of ansatz (2). 

The resuhs of the fits for SU{2) and SU{3) 
are shown in Tab.l. Although we do not expect 
a second order transition for 4-flavour, Nc = 2 
QCD in the continuum limit, we find no indica- 
tion for a discontinuity, i.e. a first order phase 
transition. This may indicate that the transition 
is still controlled by the U{1) symmetry of the 
staggered fermion action, i.e. we are still at too 



strong coupling. 



5. Conclusions 



4. Finite Temperature Potentials 

We have calculated finite temperature poten- 
tials on 16^ X 4 lattices for m = 0.025 and T < 
using Polyakov loop correlations up to a distance 
of RT = 2.0 



PLC{R) = (Ll'(0,0,0)i(dxi,dx2,dx3)} 



(3) 



with dxi = 0..4 and on-axis dxi — 0..8. 

The potentials are extracted from the following 
relation: 



(4) 
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Figure 5. Comparison of the Potentials V{R) for 
= 4 - SU{2) and Nf = 2- SU{3) for T « Tc 
and smaller temperatures. 

We compare our potentials with results for two 
flavour QCD |4|. For both theories we see clear 
signs of string breaking. The potentials seem to 
rise up to a critical distance Rc where the string 
breaks and V{R > Rc) stays constant. It seems 
that for SU{2) the string breaks at smaller dis- 
tances than for SU{2>) which may be attributed 
to the larger number of hght states in the SU{2)- 
theory that can contribute to the string breaking. 
It is however not yet clear how the heights of the 
potentials are controlled by self-energy contribu- 
tions. 



So far we have no indications for a first order 
phase transition for QCD with four flavours and 
colour group SU{2). This may be a strong cou- 
pling artifact due to the incomplete chiral sym- 
metry of the staggered action. Further investiga- 
tions with larger temporal extent are needed to 
get closer to the continuum limit. 

A rather striking feature of the SU{2) theory is 
the decoupling of chiral symmetry restoring and 
deconfining aspects of the transition. While the 
chiral order parameter changes more rapidly close 
to Tc with decreasing quark mass the Polyakov 
loop expectation values become flatter in this re- 
gion. 

Our data are compatible with earlier results by 
Kogut et al. Simulations at smaller quark 

masses should be performed to see if a first order 
phase transition appears as suggested by Kogut 

i- . 

The finite temperature potentials show signs of 
string breaking for SU{2) in agreement with the 
SU{3) data of C. DeTar et al. 0. Due to the 
larger number of light states in SU{2) it seems 
that the string breaks earlier and the plateaus 
are lower than for SU{i). 

So far this is a qualitative comparison. One has 
to analyze the self-energy contributions to the po- 
tentials and needs information about the lattice 
spacing to really compare both. 
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